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I. INTRODUCTION AND DISCUSSION OF RESULTS 

Advances in semiconductor technology have renewed interest in the properties of one-dimensional (ID) 
electron systems. It is well known that the electron-electron (e-e) interaction alters properties of ID 
system qualitatively. To gain a better understanding of the problem one should try to clarify the physical 
question of a primary importance: what is the nature of the ground state of this system? In order to get 
the answer one usually studies the ground state for ID interacting fermions using the ,, density-density" 
correlation functions. However this information is not direct. As it has been shown the correlation 
functions of the problem contain the terms that decay very slowly. Usually they are interpreted in 
the following way: the oscillations with the Fermi momentum (pf) doubled are related to the Peierls 
instability (Refsi^) while the oscillations behaving as App are interpreted as a marginal Wigner crystal^. 

In Refi^ the wave function of the ground state of spinless fermions was constructed for the exactly 
soluble Tomonaga - Luttinger model. It has been shown that at sufficiently low temperatures the system 
should be in the state that has nothing in common with a system undergoing the Peierls transition. It is 
a state with a spontaneously broken chiral symmetry. There is a long-range order in the electron system. 
In the limit of infinitely strong interaction, at low temperatures a condensate of finite density is formed. 
It consists of neutral (exciton-like) pairs of a right electron and a left hole or vice versa. The uniqueness 
of a one dimensional (ID) system impels one to consider a phase transition of the 2nd kind in a channel 
of a finite length Lu. The point is that the temperature of the phase transition vanishes as i/L^ as it 
should be. Note that actually the phase transition temperature need not be too small. For the length 
Lit ~ cm it should be about 1 K. In other words, consideration of the limit L\\ — > oo cannot be 

sufficient for description of the modern experiment. 

Traditionally many-component fermions in the ID systems have been extensively discussed in the 
literature. The interest was going to rouse by separation the spatial and spin degrees of freedom 1 , 2 . In 
the present paper we are going to discuss the form of the ground state wave function for this case. As 
the variables are separable one can expect that the ground state wave function is a direct product of two 
factors where one of the factors, describing the spinless component, coincides with the ground state wave 
function of spinless fermions. In the present paper we show that in fact one has an entirely different 
situation. Namely, the most correlated state is the state where the total spin vanishes. However, for 
n-component electrons this state consists not of pairs (as in the case of the spinless fermions) but of 
pointlike neutral complexes containing 2n particles and having a chirality ±2n. 

For the ordinary electron system, n = 2, one has the complexes consisting of two right electrons and 
two left holes with opposite spins. For the conducting nanotubes n = 4 (Refs000), and the complexes 
consists of eight particles. The complexes with a smaller number of particles can have a nonzero spin 
but their correlation is much weaker. For example, for n = 2 the spin phase can be realized only as 
a Kosterlitz-Thouless one, and for the limit of infinitely strong interaction the density tends to zero as 
1/ \/L\\. By contrast, the spinless phase has for this limit a finite density. 
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This situation is typical for many field-theoretical models with Adler anomaly (this is the case for Lut- 
tinger model as well). It is known that in such models the new fermion interaction ("t'Hooft interaction"—) 
can appear with a vortex which is a product over all components of fermions. In many cases t'Hooft 
interaction leads to the spontaneous breakdown of chiral invariance with 2n fermion order parameter. In 
particular, this is a case for multicomponent Schwinger modelS and, as we shall see, for the Luttinger 
model in the limit of infinitely strong interaction. 

For last model the most correlated state is built out of complexes each of them having the maximal 
possible number of degrees of freedom. This state differs qualitatively from the result of a common 
treatment, i. e. marginal Wigner crystal (— ). Instead of the phase transition of almost first kind one gets 
the phase transition of almost second kind. In order to manifest breakdown of the chiral symmetry in the 
Luttinger liquid we have calculated exactly the wave function of the ground state in this model, Eq. (|27p. 
and demonstrated explicitly that its symmetry is less than the original symmetry of the Hamiltonian. 
(This is the definition of spontaneous symmetry breaking) . These are basically the main results of the 
paper. 

For one component fermion system pointlike complexes with more than 2 particles are forbidden by 
Pauli principle. This is not the case for the multicomponent fermion. For this reason, expression for the 
ground state wave function Eq. (|27|l is much more cumbersome and the calculations are more involved. 
Therefore we will present the results only for a short range potential in the limit of infinitely strong 
interaction. At the same time, in our case the physical picture is quite similar to that of the one 
component fermions. For instance, had one taken into account corrections in the reciprocal strength of 
interaction one would have gotten a Kosterlitz-Thouless phase for spinless complexes too. The strength of 
the interaction going up, it would be gradually transformed into a state of a definite condensate density. 

The paper is organized as follows. In Section [fl we have given a brief review of our results. Section ITT1 
contains discussion of the main difference between the multicomponent problem and spinless one in regard 
of the theoretical description and most essential steps of calculation. In Section IlIII we give arguments 
concerning the possibility to apply the theory to nanotubes. Appendix is devoted to a derivation of some 
intermediate results. 



II. DESCRIPTION OF APPROACH AND DERIVATION OF THE MAIN RESULTS 

Our starting point will be the usual Tomonaga — Luttinger Hamiltonian H (see, for instance, Refslld 
2) for a system of interacting electrons without regard of the backscattering. For the case where the 
interaction does not change the electron spin the Hamiltonian H can be expressed through the density 
of the right (R) and left (L) electrons 

Q {x) a = Q B ,a (x) + QL.ol (x) 

(the spin index a equals ± for the spin ±1/2 respectively) as 

H = J2 dx *R,a ( x ) v f {-id x ) (x) + & L a (x) v F id x i> L . a (x) 

a J 

+ dxdyg(x)V(x-y)g(y). (1) 



Here v F is the Fermi velocity and 

* a (x) = exp (ip F x) ^ R:a (x) + exp (-ip F x) * iiCe (x) , (2) 

while g is the total density. For simplicity, we assume that the interaction is spin-independent as well (it 
is not too essential for calculations). Now we introduce the electron and hole operators in the usual way 



*(R,i)a (x) = J ^ (exp (±ipx) d {R , L)a (p) + exp fcipx) ^ {R L)a (p)j = (3) 
o 

a {R ,L )a (x) + V {R V)a (x) 
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Here a_R,+ (x) is the operator of annihilation of an electron with spin +1/2 while 6_r,+ (x) is the operator 
of annihilation of a hole with spin —1/2. To write such interaction via functional integral it is necessary 
to introduce one Bose field $ and apply a version of the Hubbard — Stratonovich identitj*ii: 



exp 



X>$exp 



■' X 

T 



■ T roc i 

dt I -£v(p)g{ P ,t)g(-p,t) 

2 Jo J-oo 27T 



' dt I ^ (e (p. *)*(-p.*) + e(-p. *)*(?.*)) 



Here V (p) is the Fourier-transform of the e-e interaction and the normalization factor is equal to 



A/"= /x>$exp % - ( dt [ ^-^(p,t)^(-p,t)V- 1 
J 2 Jo J-oo 2lT 



(P) 



(4) 



(5) 



This identity shows that the theory with e-e interaction is equivalent to a theory of non-interacting 
electrons in an external field (Had one taken a theory with spin-dependent interaction he has to 
introduce two independent Bose fields. The calculation would be more cumbersome but the physical 
picture would be the same.) 

Calculation of the ground state wave function is given in detail in Ref^. It is based on calculation of 
the evolution operator for the electrons 



m >< m\. 



(6) 



Here \n > are the exact wave functions of the Hamiltonian H in the secondary quantization representation, 
T is the observation time. S(T) determines the evolution of an arbitrary initial wave function (< m\) 
from the time t = up to final states \n > (at t = T). (Henceforth we imply that the Schrodinger 
representation for operator with time-dependent wave-functions is used.) 

Eq. © suggests the general method to obtain the wave functions. One has to calculate first the 
evolution operator and present it as a sum of time-dependent exponents. The coefficients in front of 
these exponents are products of the exact wave functions and their complex conjugates. In order to 
extract the ground state wave function one has to take the limit T — > oo (we add infinitesimal imaginary 
part to the energy). Proceeding to Euclidean time (T — > —i/Q) we see that evolution operator determines 
the density matrix for the equilibrium system at a non-zero temperature 0. 

As for the case of spinless electrons (see^ for the details), one can present the evolution operator for 
the electrons in an external field as 



S ($) = exp (So + log[Det $(T)])\F><F\ 



(J) 



Before integration over the fields $ the equation for S ($) undergoes some changes in comparison with 
the spinless case. They amount to appearance of a factor n, the number of components of the electron 
wave function, in the equation describing the quantum fluctuations in the electron system under the 
action of the field $ (T) (in the spin case n = 2) 



In [Det $ (T)) 



n 

ill 



dtdti 



dp 
2^ 



$ (-p, t) $ [p, t x ) \p\ exp [-i\p\v F \t - t x 



(8) 



The operator structure of Eq. J7J) is determined by the second part of action, S$. Here one should take 
into account the spin index 

So= ! dxdx ' \ka (*') G " 0; x, e) a i<a (x) + a\ a (x') G° (x', T;x,T- e) b\ a (x) 



i—R.L:a 



d\ a (x') G? (a/, T; x, 0) a i>a (x) - b i<a (x') G° (x\ 0; x, T) b\ a (x)} , (9) 
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where Gf is the spin-independent free particle Green function 

G° R L (x, t; x u ti) = [v f (t - ti) T (x - x x ) - iS sign {t - ti)]' 1 . (10) 

ZIYI 

Behaviour of the multicomponcnt fermions in the external field is quite similar to the one-component 
case. Essential complication, however, appears after the integration of the operator S (<E>) over <E> with 
the weight 

Two points should be indicated. 

• The coefficient n in Eq. © that enters the equation for the action. Due to this coefficient, af- 
ter calculation of the integral over $ the analytical structure of the resulting expression becomes 
much more complicated. As a result, in the integrals defining wave functions of the multiparticle 
complexes, one gets cuts instead of simple poles, as in one component case. This leads to rather 
cumbersome complex wave functions. In particular, complexes are non-local. 

• The nonlocal property brings about essential enhancement of the various electron states. A number 
of electron states is forbidden for one component electrons due to Pauli principle. In contrast to 
this, in the multicomponcnt case the number of connected diagrams becomes infinite. This makes 
the expression for the ground state wave function rather cumbersome. However the very fact of 
existence of chiral phase for the infinitely strong interaction still persists. 

Because of non-Gaussian form of the final functional integral it is impossible to perform the integration 
in Eq. (JSJ over &(x,t) in a closed form but it is possible to obtain an arbitrary term of the evolution 
operator expanding it in Sq . It will be sufficient in order to get the ground state wave function as we will 
have the integral of a Gaussian type that can be easily calculated. After doing the final integration over 
$ the final recipe of calculation of the evolution operator can be written as a sum of the following terms 

So (W), ■■•■) ex p ( s f )\f >< n (ii) 

Here \F > is the Fermi "sphere" while the terms in Sq determine the operator structure of the wave 
functions, i. e. all the possible particle configurations as a result of their interaction. (|ll|l is a sort 
of symbolic expression. Indeed, the analytic equation for the effective action S° s in the nth term of 
expansion depends explicitly on the particle configuration in the preexponential factor (<So) . Naturally, 
it is different for different terms. Note that the evolution operator 10 is determined in such a manner 
that the initial state expressed through the electron and hole annihilation operators and the final state 
determined through creation operators are given at different times. It means that during calculating the 
evolution operator one should consider the operators a R L (x) and an,L (y) as anticommuting. 

In order to write the expression for 6>° ff , we introduce the following notation for the coordinates of the 
electron-hole creation-annihilation operators: 

1. We will denote by x (y) the coordinates of the right (left) particles. 

2. We will put a tilde on the coordinates related to annihilation operators (the initial state): the 
coordinates of creation operators (the final state) will have no tilde. 

3. We will prime the hole coordinates. 

The effective action differs from the action for the one-component fermions only by a factor and in the 
limit of strong interaction 

m y>1 

ttvf 

is equal 

s cS = — - 2J I — | [Rf (-p, xi ■ ■ .) fcf 0, xx . . .) + Tii (-p, xi . . .) (j>, x-l . . .)] 
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7En ex P ( J #^) n f (-p, *i • • ^ Cp, «i • • •) , (12) 



The extra factor n is the number of the fermion components. It appears because in the Luttinger 
model the excitation spectrum isi£: 



Up = \p\ VF Jl + ?¥M. (13) 
y 7tvf 

Eq. (|12[) is valid in the temperature region 

©chiral « 6 « C = w (p»m)- 

The right-hand side of the last inequality is the energy of excitations with a minimal momentum (for 
periodic boundary conditions p mm — 2w/Ln) while S c hiral = |Pmin|vF is the degeneracy temperature of 
the ground state. In this temperature region one can neglect the energy difference between the states of 
different chiralities. This means that the ground states with different chiralities becomes degenerate. 

The origin of this inequality has been discussed in detail in Ref. 4 . It is not sensitive to the number 
of the wave function components. For the temperatures <C Ochirai the last term in Eq. (|12H should 
be omitted. Then the corresponding equation will be valid for the low temperatures too (see detailed 
discussion in Ref. — ) . 

The functions TZij (p, x\ . . .) in the equation l|12f> depend explicitly on the electron (x..) and hole (x...) 
coordinates in the preexponential factor. These functions are given by 

Kf(p,x...)= ^2 (p) [exp (ipx) 

x...;x'...;y...;y'... 

- exp (ipx')] + 9 (-p) [exp (ipy) - exp (ipy 1 )] , 
TZi(p,x...)= 9 (-p) [exp (ipx) (14) 

x...;& ... 

— exp (ipx')] + 9 (p) [exp (ipy) — exp (ipy 1 )] 

To obtain the complete expression for the ground state wave function one has to consider all the 
complexes and separate the connected parts out of them. It is not necessary, however, since, according 
to the general theorem, the complete wave function is the exponent of the connected complexes^. 

For one-component fermions there is only one possible 2-particle connected complex. This is not true 
for multi-component case, many of the scattering channels are possible, so the number of connected 
diagrams is infinite. In principle, one can calculate the exact wave function of any given complex taking 
the Gaussian integral in $. Unfortunately, it is not enough to present the whole wave function of the 
system in a closed form, but actually in order to prove the existence of symmetry breaking we do not 
needed it. To be sure of the fact it is sufficient to prove that the wave function symmetry is less than 
Hamiltonian. In order to do this one should analyze the simplest connected diagrams bringing about a 
spontaneous breaking of the Hamiltonian symmetry. The rest terms either have the symmetry of the 
Hamiltonian or describe the scattering of the simplest correlated complexes and also violate the chiral 
symmetry. 

Now we embark on analysis of the simplest diagrams of the evolution operator for the electrons having 
a spin. We will begin from the temperature region c hirai 3> @- In this case the action and therefore 
the evolution operator factorize so that one can consider explicitly the ground state wave function |fi >. 
The simplest nontrivial diagram we are interested in is 

dx+dx' + dx-dx'_ dy+dy' + dy_dy'_ d R , + (x+) b R + (x' + ) a R _ (x-) b R _ (x'_) 
(2mf (2mf x' + - x+ - id x'_ -X--i6 

dl + (y + )bl + (y' + ) a* ,_(!,_)&£,_ (y'_) 



y+ -y+-iS y- - y' 



exp[5 / ff ( 2 ; + ,...)]|F>. (15) 



Further on we will see that the terms with a smaller number of operators give a weaker correlation than 
Eq. dSJ. 
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The action for this electrons-holes configuration is 



1 U a , a >... (Sa ~ V/3 + **) 6v ~ Vp ■ , 

S f eii (x, a ,...) = -ln ■ ) f. (16) 

" U a ,a'... Oa -V0 + (^a' - y'f}> 



id 



In fact, it differs from the corresponding expression for one-component fermions by the factor 1/2. This 
factor leads to the singularities of the integrand Eq. (|15|) which are cuts instead of simple poles in one 
component case. This prevents us from explicit calculation of the integral. 

Nevertheless, it is possible to recognize the spontaneous breakdown of the chiral symmetry in our 
system. To do this, several steps are necessary. First, we have to analyze what new bound complexes 
appeared as a result of interaction. We have to take the arbitrary connected diagram and try to separate 
complexes with a smaller number of particles out of it. To do this one has to consider all particles in one 
complex as being one close to another, while the distances between different complexes are large. If the 
full wave function in this limit decays into the product of two wave functions, one depending only on the 
coordinates of first complex, the other only on coordinates of the second one, then the complexes can be 
considered as new " free particles"—, as the probability to find one such complex does not depend on the 
position of another one. In other words, we should present a term of expansion of evolution operator we 
consider as a product of the form: 

dx + dy + ± ± 

■ ■ ■ ■ a k+ i x +) ■■• K ( x +, ■ ■ ■ > v+ ■ ■ ■ ; ) (y+) ■■■ ■ 

Then we should check that provided the complex (x+, . . . ,) is moved off from the complex (y + . . . ;) over 
the distance of the order of L\\ the amplitude K (x + , . . . ,y + , . . .) tends not to zero (as is usually the 
case with scattering amplitudes) but to the factorized product k (x+, . . . , ) k\ (y+, . . .). Here each factor 
depends on the variables of the first or second group. This means that two complexes are formed as 
a result of the interaction. If the intercomplex distance is large enough, their contribution to the wave 
function can be presented as follows 

.k(x + ...,) at + (x+ ...) [ -^=fci (y+ ■ ■ ■ , ) fit + (y+ ■ ■ •) \F > . 



The remnant part of K (which is K — k ■ ki) is a connected diagram which describes the intercomplex 
scattering. The theorem of logarithm connectedness 13 guarantees that the same connected complexes 
will appear in the next orders as well with correct combinatorial coefficients and the final answer is an 
exponent of connected complexes. In particular, the first order term 

■ k(x + ,...)a f R + (x+)...\F> 

should appear in the expansion of evolution operator directly, unless it is forbidden by some conservation 
law (e.g. chirality conservation for the lowest temperature case in our model). In this case one has to 
use the projector on the proper state as in (Eq[^J). 

Whether a phase transition of the second kind exists or not, depends on the symmetry of the complexes 
ajj , (x + ) ... If they are less symmetric than the initial Hamiltonian one has a symmetry breaking. As 
a result it is possible to introduce a non-vanishing order parameter in the less symmetric phase while 
it vanishes in the more symmetric one where the complexes do not exist. (More precisely, taking into 
account a fluctuations of the low symmetry phase in phase with non-broken symmetry one can state that 
an order parameter should not enhance with Lu in the last case. Our definition the order parameter is 
given in Eq. I|29(1 .'1 According to Landau (see Landau and Lifshitz l i) this is the definition of the phase 
transition of the second kind. If, however, the symmetries of all the connected complexes and of the 
Hamiltonian are the same one has a long-range correlations without spontaneous symmetry breaking. 

Now we embark on analysis of the simplest diagrams of the evolution operator. By analogy with the 
theory of one component fermions one could assume that the simplest connected diagrams originate from 
the term 



r dx + dx' + dy_dy'_ (x' + ) (g-) b\ _ (y'_) ^ 

J 2ni 2ni x' + — x+ — iS y- — y'_ — id cff 
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However, because of the factor n the corresponding contribution to the action is 

5/ (a; , , .) _ 1 i n (g± ~ ^ + ^ « ~ ^- + ^) , (is) 
c ' 2 (a;^_ — y_ + ifi) (x + — y'_ + iS\ 

The bound chiral complexes are determined by the singularities of the integrand at \x' + — y_| ~ 
d, \x+ — y'_ \ ~ d, \x + — y_| ~ L11 (here d is the width of the conductor). As a result, the contribution we 
are interested in is of the order 

dx + dy-a) R (x+) b[ (x+) a[ (y_) b R : - -\F > . 

\ x + V—\ 

This quantity tends to at \x + — — > 00, but more slowly than for free particle^ case. 

Now we will show that the most correlated state can be obtained from Eq. ( 115(1 . It can be split into 
two four-particle complexes, each with zero spin, having the chiral charges ±4 (a R + a R _b L V* L and 

aJ L + a) L b R + b R _). The amplitude K in Ea.f I15fl factorizes and does not tend to at — > 00. Indeed, 
the c-factor in the integrand is 

1111 

K(x + , ...) = - 



x+ — iS x'_ — x- — iS y+ — y' + — id y_ — y'_ — id 



Ua,a'... ( x a ~ Vli + £<0 [x' a , ~ y' p , + i6 

)jn a , a >... « - ^8 + ^) (x a , - y' fj , + i5^j 



(19) 



Now, let us consider the regions of integration x + ~ X- ~ y'_ ~ y' + and x' + ~ x'_ ~ y+ ~ y_ assuming 
that the distances between these groups of variables are of the order of Lu. Then the amplitude K tends 
to 

V +A (x+,...)V- 4 «,...) = 



= 1/^ (x + -y' + + iS) (x + -y'_+ iS) {xi - y' + + id) {xl - y'_ + iS) 

(x' + -y+ + iS) (x' + -y-+ iS) (x'_ - y+ + iS) (x'_ - y_ + iS) , (20) 

This means that each amplitude V depends on the variables belonging either to the first or to the second 
group. This property of the amplitude permits one to single out of the full equation for the evolution 
operator the connected complexes and the amplitude of intercomplex scattering that tend to for large 
intercomplex distances. 

This term of expansion, besides the chiral complexes, has also neutral ones with zero chirality, 
®R J*L _ ana - ®r -b R _a\ b\ ■ They do not violate the symmetry of the Hamiltonian. How- 

ever, they should be singled out, so that the remaining scattering amplitude tended to zero in the whole 
region of the variable variation. This permits one to interpret it as the intercomplex scattering ampli- 
tude. The zero chirality complexes are not important for existence of the phase transition. However, 
one should take them into account for calculation of the matrix elements as the theory has no small 
parameter to neglect them. To check that they exist we consider in (|15f) the region x+ ~ x' + ~ y- ~ y'_ 
x- ~ x'_ ~ y+ ~ y' + . In this region the amplitude K tends to Vq (x+, . . .) Vo (x-, . . .) where 



1 \ (x a - y- a + id) (x' a - y'_ a + iS) 

Vo (x a ,...) = — } -r V (21) 

(x' a -x a - iS) (y_« - y_ a ~ iS) _ yl _ a ~ 7q ^ _ yl _ a ~ ^ 

This means that this quantity can be presented as a product of the amplitudes, each of them remaining 
finite provided the distance between them tends to infinity. 

Now it is convenient to introduce the amplitude of intercomplex scattering V co u. By derivation, it 
tends to zero, the intercomplex distance tending to 00: 

VcM (x+,...)=K(x+,...)- V +4 (x+, . . .) V- 4 (x' + , ...)-V (x+,...) Vo (x., . . .) . (22) 
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The contribution to the ground state wave function we are interested in can be presented through these 
amplitudes as 

dx + dx' + dx-dx'_ dy + dy' + dy-dy'_ , - f t 

a R,+ ( x +) b R,+ ( x +) a R,- ( x -> 

{2m) {2m) 

&k- 00 < + <y+) 04) a L (»-) *L (y-) 

(V +i {x+, . . .) y_ 4 «, • • ■) + V (x+, . . .) V (!_,...) + V co « (x+, . . .)) |F > • (23) 

First term here describes the non-interacting complexes with non-zero chirality, second describes the 
4-particle neutral complexes and third (connected part) is related to their interaction. We are mainly 
interested in first term as it is connected with breakdown of chiral symmetry. 

The chiral complex which we obtained is already a connected one and cannot be separated to simpler 
ones. This means that its wave function is a decreasing function of interparticle distances. It is shown in 
the Appendix that probability to find particles of the complex far from each other is negligibly small. 

Besides, one should take into consideration that in the temperature interval 8 <C Ocimai where one 
need not consider the last term in Eq. (|12|) there is one to one correspondence between the complexes 
with chirality Q = 4 and Q — —4, so that the total chirality of the state equals zero. The theorem of 
logarithm connectedness states that the ground state wave function can be presented as 



1° >= Pc^oexpy^Tr 1 V+ 4 (x a , £_ a , y'_ a , y' a ) (a R a {x a ) d R _ Q {x_ a ) 
{2m) v 

X %.,-a (y'- a ) bl, a {y'a) + 4,a (t/J (v'-o) Vr,-c (»-«) &L (x a )) 

+ /n 1 .,2 V +0 (xa,X-a,y'- a ,y a ) (a Ra {x a )b R (x' a ) 

{2m) v 

X «L,-a (!/-«) bl,- a (y'-«) + 4,-a (*-«) ^R.-a W «L &£,a 

+ , 1 .,4 ^°» (»a, X-a, ■■■) (la) &^ (x_ a ) ff L (y'_ a ) P L a {lj' a ) 

{2m) 

X 4,a (»a) 4,-a &k,-a Ka) OO + " ' '] ' ( 24 ) 

Here Pc=o is the projector on the state with zero chirality. The symbol Tr includes the integrations in 
over particle coordinates. The terms omitted in Eq. 124fl describe scattering of three and more complexes 
while all the elementary complexes are present here. Mark that the complexes with non-zero chirality 
have appeared in the theory. Nevertheless, the wave function of the ground state as a whole described the 
state with Q = in this temperature region, i. e. the symmetry of the ground state is the same as that 
of the Hamiltonian. The states with nonzero chirality have a bigger energy (of the order of 2ttvf/Li\). 
Therefore the spontaneous symmetry breaking may take place only in the region of higher temperatures 
6 6 c hirai where such an energy difference is not essential. In this temperature region one should 
consider in the equation @ for the action Sq also the term where the time arguments of the Green 
functions differ by T. (Practically it is more convenient to introduce the temperature by replacement 
T — > — i/0 in the final equations.) Then one has the following non-trivial term in the operator of evolution 

1 at + (x + ) a R (x +) at_(x-)a R (x-) 

Ir- 



{2m) x + — x + + VfT — id aL — + VfT — i5 

— exp[S^ s ]{x + ,...)\F >< F\. (25) 



V+-y' + - VfT + iS y'_ - y'_ - v f T + id 



The action for this configuration is 



i U a<a >... (y'a -v'p- v f T + iS ) (a?o' - x f}' - v f T + iS ) 

In ^ <- . (26) 

n. a ,a>... \ X <* -y'p + l5 ) (y'a' - x 0' + l5 ) 
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It is readily seen out of the operator structure of this term that the amplitude V4 appears here au- 
tomatically (without extracting the amplitudes of the neutral complexes and channels of scattering). It 
is a consequence of the theorem of logarithm connectedness. It guarantees coincidence of the amplitude 
in this term with V4. Indeed, in the region where the same variables with spin up and spin down are 
quite near to each other (for instance, x a ~ X— a ) while the coordinates in the creation and annihilation 
operators are apart at the distance of the order of Lu (x a ~ x a ~ L 3> 2"KVf/&) all the c-factor in the 
integrand of H25|) turns into the factorized expression 

V+4 (x+, . . .) F* 4 (x+, ...). 

This proves that the chiral four particle complexes were singled properly out of a more complicated 
expression (|22|) . Such terms in the evolution operator result in any chirality of the ground state. However, 
any state with a fixed chirality should be unstable relative to the backscattering, however weak, violating 
the chirality. Therefore one should consider a superposition of all such states. As a result, one has, in 
the same way as in the theory of superconductivity, to introduce a condensate with a fixed phase instead 
of a state with a fixed chirality: 



|fi > e = expJ^Tr 



2^+4 (i«> x-a, y'- a , y' a ) (exp (iO)a Ra (x a ) a R _ a (x_ Q ) 



(2m 

X bl- a (y'-a) b f L>a (y'a) + exp H0)4,a (Va) (v'-a) ^R,-a (*-a) (*„)) 

+ 7^2^+0 (xa,%-a,y'- a ,y' a ) (a R , a (x a )b R (x' a ) 
(2m) v 

X 4,-a (W-o) tt L ,-a (y'-a) + (*-«) ^R,-a ( X '-a) A,a (»«) &i,a (fa)) 

+ 7^4 ^ -■) a k-« ^i,-a (V'-a) !>L. a (Va) 

(27Tl) 

X al,a (Va) «L- Q (l/-a) S^_ a (^) + ....] , (27) 

Eq. (|27(l demonstrates that for a strong electron-electron interaction a spontaneous symmetry breaking 
takes place. Two first terms in the equation for the ground state wave function are not invariant under 
the chiral transformation 

®r,l (x) -» exp (±iA)y RjL , (28) 

while the Hamiltonian still retains the invariance. (Here A is a constant.) 

Such a form of the bound state permits to introduce the order parameter. It equals to zero (or, to be 
more exact, not grows with Ln) in the phase of high symmetry (O 3> C = w Pmln ) and is proportional to 
Lii in the phase of low symmetry for low temperatures. (One can prove the first statement using ordinary 
symmetry considerations or by direct analytical calculation.) For our case the following quantity can be 
considered as the order parameter 

dx e < n|4_ Q (x) a R _ a (x) b\ a (x) b\ _ a (x) \Q > e ~ Ly . (29) 

One can see from Eq. I|29(l as to why the phase transition of the second kind demands the chirality 
degeneracy of the ground state. For the order parameter to be nonvanishing it is necessary to be able 
to add to the ground state an extra four particle complex. This demands condition O 3> ©chiral because 
zero chirality state has the lowest energy and the energy difference between the state and other ones 
about ©chiral- There is also the upper bound of existence of the low symmetry phase 

©<w Pmin . (30) 



This limitation is due to the fact that in one dimensional systems the long range order is suppressed by the 
thermal excitations. In the temperature region given by Eq. (|30|l one can neglect it; this condition does 
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not depend on the number of the fermion components and is discussed in detail in*. So, the temperature 
region 

Ochiral < 6 < 9 C (31) 

is the region where the broken phase exists. 

Thus in the system of interacting multicomponent fermions the most correlated state consists of 2n 
operators and has the chirality ±2n. It is this state that in the limit of infinitely strong interaction results 
in the phase transition of the second kind. At the same time the spin complexes of a smaller number of 
operators, violating the symmetry of the ground state, can exist only as a Kosterlitz — Thouless phase. 
For the limit of infinitely strong interaction their density tends to zero as By contrast, the 

spinless phase has for this case a finite density. This rule is quite general. 



III. CARBON NANOTUBES 



The conducting carbon nanotubes give one more example of multicomponent electrons. In order to 
generalize the theory developed above one should have one dimensional conducting tubes and such e-e 
interaction that could be rewritten in the density-density form — see Eq.Q. It means that one should 
be able to neglect a backward and inter-component (see below) scattering. Following Ref. ^3 one can 
visualize a nanotube as a cylinder constructed of a monoatomic layer of graphite. The latter has a lattice 
of adjoining regular hexagons, so that the angle between the neighboring basis vectors, n a and rib is 
27r/3. Choosing the coordinates £i and £2 in such a way that axis 0£i is parallel to a while axis 0^2 is 
perpendicular to a one can present these vectors as 

a = a(l,0), b = a(-l/2,v / 3/2). (32) 

a is the lattice constant that is equal to <i\/3, d =1.44 A being the interatomic distance^. 
The circumferential vector L can be written as 



L = n a & + ribb. (33) 



Here n a and rib are integers. 

The electron effective Hamiltonian for a graphite sheet is 

It can be expanded in the vicinity of the points 

P = (47r/3a)(-l,0), P' = (47r/3a)(l, 0) (35) 

up to the first power in the small deviations p and p' from the values given by first and second Eq. Ij35(l 
respectively 

h(P, p) = 7 e^ (p ± - ip z ) , h(P', p') = 7 e <9 (~p' ± - r P ' z ) ■ (36) 

Here 7 = (-\/3/2)7oa; <y ~ 3 e V (see Il7ll8j) is the transfer integral between the neighboring it orbitals 
while 9 is the angle between vectors L and a. The subscripts z and _L refer to the components of p 
relative to the direction of L, namely p z _L L and px||L, so that p z (p±) is parallel (perpendicular) to the 
axis of the tube. 

The spectrum near the point P is given by 



£(P,p)=± 7v /^+Pi (37) 

where the upper (lower) sign corresponds to the conduction (valence) band in this equation. The spectrum 
near the P' point is obtained by the replacement p — * p'. 
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The electron wave function \P(r) should satisfy the cyclic boundary condition 

*(r) = *(r + L), (38) 
so that the discrete values of p± and p' ± are given by (see Ref. Il7h 




Here n = 0, ±1, ±2, . . . while |L| = aW n\ + — n a n,b m , v = or ±1 and is determined by the presentation 
of the sum n a + rib as 3N + v (N being an integer). The nanotubes are conductive (metallic) for the 
combination 

n = v = Q (40) 

and we will consider this case for our further analysis. In other words, in such tubes there are two conic 
bands, i. e. the points apP with ap = ±1. The big phase corresponding to the momentum apP should 
be extracted in the same way as it has been done in Eq. Besides, we assume that due to the presence 
of gate electrodes the Fermi level is well above (or below) the points apP (cf. with Ref£). As a result, 
we will have a theory with four-component fermions. There are two extra branches corresponding to two 
values of ap. In each of them there are analogues of R- and L-particles. To introduce them one should, in 
full analogy with Eq. (J2J, separate the corresponding phase factors with large phases. In order to establish 
correspondence between the present model and the Luttinger one we should be able to neglect both the 
transitions between different branches (different values of ap) and within the same branch between R- 
and L-particles. As indicated in Ref. lltA the nanotubes have comparatively large radii that encompass 
with (N 3> 1) atoms. Therefore the only e-e scattering that is important in this limit is the forward one 
that involves a small quasimomcntum transfer. The matrix element describing the backscattering within 
a band as well as the P «-> P' scattering acquire an extra small factor of the order of 1/N. This is why 
one can neglect these types of scattering. This means that we can use the results obtained in Section ITT1 
Repeating the arguments of this section for n = 4 one can come to the conclusion that a condensate is 
formed in the ground state. It consists of the eights of the form 

& R, a , ap fa) *R,-a,a P (») ^,a,ap (*) b l,-a,a P ( x ) & R,a,-a P ( X ) & R,-a,-a P ( X ) &£,a,-ap ( X ) ^L,-a,-a P ( X ) » 

Their chirality is ±8. 
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APPENDIX A: CHARACTERISTIC DIMENSIONS OF CORRELATED COMPLEXES 

The simplest way to give the proof of existence of a bound chiral complex and to determine its char- 
acteristic dimensions is to consider the state with a single chiral complex: 



|* c >= Tr 



4 + (z+)4-(z-)£L-(y-)SL + (2/;) 



n 



a.,a'=± 



(x a - y' a , + i5){x a - y' a , + iS) 



\F>. 



(Al) 



For instance, one can calculate the probability to find an electron with spin up at the distance \z+ — z_ | > d 
from the electron with spin down as 



< 



$ c |&fl,+ (z+) a R .+ (z+) a R _ (z~) a R ,- (z_) |$ c > . 
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Moving all the creation operators to the right and all the annihilation operators to the left one gets 
A 2 /\z + — z-\ 2 . Here A is a constant equal to 



(o 



dx (1 



-1/2 



{l-xf + 6 : 



-1/2 



This means that the most probable is the particle configuration where \z + — z_ | ~ d (to get the physical 
parameter one should replace, as usual, 5 by d). That is the right and left electrons with the opposite 
spins are in fact always near one another forming a chiral complex. In this sense the chiral fours are 
point-like entities as the i?L-pairs for one component fermions. In the same way one can give estimates 
of the distances between all the particles belonging to a four-particle complex. 

The same calculation for a neutral four-particle complex arc a little more cumbersome. The state with 
a single neutral pair is described by the wave function 



l$o >=Tr 



(x-)b R _ (x'_)d[ + (y^ 



(x'_ 



i$) 1 (ZM 



y'+ 



isy 



(x- - y + + iS) (x'_ - y' + + i5) 



i,+ [y+) 

F > 



— y' + + i<5) (x'_ — y+ + i<5) 
In order to find the characteristic size of a neutral complex consider the matrix element 



(A2) 



It is equal to 



Tr' 



< $o\a/ R _ (z-)a R - (z_) d\ + (z + ) a L+ (z + ) |$ > 



iS)~ 



(x'_-iS) l (-y' + -iSy 



y'+- y-i 


. -i5) 1 (W-+i5) 1 (- 


y'+ + 


isy 1 




z + — iS) (z_ — z + + x'- 


-¥+ 


-iS) 




z+ - y' + - iS) (z_ - z + 


+ S 7 _ 


-iS) 




z + + iS) (z_ — z + + x'_ 


-y'+- 


\-iS) 




z+-ff++ i5) (z- - z + 


+ x'_ - 


f iS) 



(A3) 



(Tr' implies that one should integrate over all the variables besides z a ). The exact expression for this 
matrix element for arbitrary values of AZ — z + — z_ is rather cumbersome and non-informative. It is 
sufficient to prove that the most probable is the particle configuration where AZ ~ d. 

To do this let us note that for AZ ~ S all the integrals in Eq. (a3) converge and are dominated by the 
regions of the order 5. Let us show that for A2 > i5 the matrix element (|A3(1 has an additional small 
factor 8/ AZ. Consider, for instance, the integration over x'_. Only the factor 



(z- — z + + x'_ 



y'+ 



iS) 



^ [z- - z + + x'_ -|- i5) 



(A4) 



in the integrand has a singularity in the lower semiplane. The rest integrand does not depend of AZ 
and has singularities only in the upper semiplane. In the main approximation in S/AZ HA4|) tends to 
1. (AZ 3> S, while the regions x'_ , y' + , giving the main contribution into the integral ~ i5). Hence in 
this approximation all the integral l|A3(l vanishes. It is nonzero only in the next approximation due to 

the factors of the sort J (x'_ + iS^j / (z_ — z + ) <C 1. Thus the probability to find the electrons we are 
interested in at a large distance is small. Most probable is a four where these particles are at the distances 
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